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Introduction 



47G20 Integro-differential operators, 34L40 Particular operators (Dirac, one-dimensional Schrodinger, etc.) 
35K90 Abstract parabolic evolution equations, 35L90 Abstract hyperbolic evolution equations. 

Introduction 

A canonical form of many linear evolutionary 1 equations of mathematical physics is given 
by a dynamic system of equations 

d V + AU = F 
completed by a so-called material law 

V = MU. 

Here do denotes derivation with respect to the time variable, A is a usually unbounded 
operator containing spatial derivatives and M. is a continuous linear operator. Here we 
would like to inspect more closely a very specific situation, where the space dimension is 
1. Although this is a rather particular case it has the advantage that an impedance type 
boundary condition, which we wish to consider, can be considered in a more "tangible" 
way without incurring regularity assumptions on coefficients and data. In the higher di- 
mensional case for the acoustic equations, which can also be discussed with no further 
regularity requirements, the constraints on the impedance type boundary condition are 
much less explicit, compare [7], [5]. Moreover, we hope to gain a different access to a class 
of problems, which are closely linked to Sturm-Liouville operators, in fact yielding a gen- 
eralization of such operators. That we are discussing the direct time-dependent problem 
rather than an associated spectral problem will actually be advantageous, since it provides 
a simpler access to the discussion of well-posedness. 

More specifically we want to consider 

/0 d\ 
A= 
\d00j 

as a differential operator on the unit interval ] — 1/2, 1/2 [ with an impedance type boundary 
condition of the form 

d a (+1/2 + 0) s ( ■ , +1/2 + 0) - v ( • , +1/2 + 0) = 
holding on the real time-line M as a constraint characterizing I w in the domain D (A) . 

w 

Here d denotes the spatial derivative and a is a coefficient operator specified more precisely 

1 We prefer the term evolutionary equations, since the term evolution equations is usually reserved for a 
special case of the class of evolutionary equations considered here. 
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later. We shall focus here on the time-translation invariant, i.e. autonomous, case. This 
means that time-translation and consequently time-differentiation commutes with a, M. 
and A. 

Our discussion is embedded into an abstract setting, which we will develop in Section 1 
first. In Section 2 we will then discuss our problem of interest as an application of the 
solution theory in the abstract setting. 



1 The Abstract Solution Framework 

1.1 Sobolev chains associated with the time-derivative 

A particular instance of the construction of Sobolev chains is the one based on the time- 
derivative do. We recall, e.g. from [4, 5], that differentiation considered in the complex 
Hilbert space H„ tQ (R) := {/ 6 L 1 2 oc (M)|(x ^> exp(-vx)f(x)) e L 2 (R)}, v e R\{0}, with 
inner product 

{f,g)^(f,g)u,o-= f(x)*g(x) exp(-2vx) dx 



can indeed be established as a normal operator, which we denote by do tU , with 

9fo d 0> „ = v. 

For 3m do jV we have as a spectral representation the Fourier-Laplace transform C y : H Ut o(M.) - 
L 2 (R) given by the unitary extension of 

Coo(K) c H„fi(R) —*■ L 2 (R) 

cf) h-> ( x h> —= I exp (-ixy) exp (-vy) (f>(y) dy) . 
V V27T J R / 

In other words, we have the unitary equivalence 

3md ,„ = C~ lr m C v , 

where m denotes the selfadjoint multiplication-by-argument operator in L 2 (R). Since 
is in the resolvent set of do iV we have that 3q v is an element of the Banach space 
L (i7„ (R), H Vt o(M.)) of continuous (left-total) linear mappings in H V) o(ti&). Denoting gen- 
erally the operator norm of the Banach space L (X, Y) by || ■ we § e ^ ^ or 



-i 



\d-H = - 



Not too surprisingly, we find for v > 

(5 »(x)= f <p{t) dt 

J— co 
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1 The Abstract Solution Framework 



and for v < 

(d »0) = - if {t) dt 



for all <p e H V fi (R) and jet. Since we are interested in the forward causal situation, we 
assume v > throughout. Moreover, in the following we shall mostly write do for do )U if 
the choice of v is clear from the context. 

Thus, we obtain a chain (H v ^ (R)) fcgZ of Hilbert spaces, wherei?^ (R) is the completion 
of the inner product space D (<9q) with norm \-\ vk given by 

Similarly, for im + v as a normal operator in L 2 (R) we construct the chain of polynomially 
weighted L 2 (R)-spaces 

with 



(R) := {/ e Lf oc (R) | (im + vf f E L 2 (R)} = H k (im + v) 



for k e Z. 



Since the unitarily equivalent operators do,u and im + z/ (via the Fourier-Laplace transform) 
can canonically be lifted to the X-valued case, X an arbitrary complex Hilbert space, we are 
lead to a corresponding chain {H v ^ (R, X)) keZ and [L\ (R, X)) keZ of X-valued generalized 
functions. The Fourier-Laplace transform can also be lifted to the X-valued case yielding 

H Vtk (R,X) ^L 2 fc (R,X) 

as a unitary mapping for k e N and by continuous extension, keeping the notation C v for 
the extension, also for k e Z. Since L v has been constructed from a spectral representation 
of 3m<9o,i,, we can utilize the corresponding operator function calculus for functions of 
3m do,!,. Noting that do = i3mdo^ + z/ is a function of 3mdo tU we can define operator- 
valued functions of do- 

Definition 1.1. Let r > ^ > and M : Bc(r,r) — > L(H,H) be bounded and analytic, 
i/ a Hilbert space. Then define 



M (do" 1 ) := £• M (-^—) C v 
v y \ im + i/ / 



where 



for (J) e Coo (R, H). 



6 



1.2 Abstract Solution Theory 



Remark 1.2. The definition of M^d^ 1 ) is largely independent of the choice of v in the sense 
that the operators for two different parameters V\,v<i coincide on the intersection of the 
respective domains. 

Simple examples are polynomials in d^ 1 with operator coefficients. A more exotic example 
of an analytic and bounded function of dy 1 is the delay operator, which itself is a special 
case of the time translation: 

Examples: Let r > 0, zv > /i e R and u e H U:0 (R, X). We define 

ThU := u( ■ + h). 

The operator th s L(H Vt o(M.,X), H u _o(M., X)) is called a time-translation operator. If h < 
the operator Ty t is also called a delay operator. In the latter case the function 

B c {r,r) b z h-> M(z) := exp(z^ 1 h) 

is analytic and uniformly bounded for every r e M >0 (considered as an L (X, X)-valued 
function). An easy computation shows for u e H u0 (W,X) that 

u( ■ + h) = £*exp((im + v)K)L v u = M(3q 1 )u = exp((^ " 1 ) 1 h) u. 
This shows that 

r h = exp (h/dQ 1 ) = exp (hd ) . 

Another class of interesting bounded analytic functions of Oq 1 are mappings produced by 
a temporal convolution with a suitable operator-valued integral kernel. 

1.2 Abstract Solution Theory 

We shall discuss equations of the form 

(d M (d^) +A)U = J. (1) 

where we shall assume that A and A* are commuting with do and non-negative in the 
Hilbert space H v q (R, H), H a given Hilbert space, in the sense that 

nt(u\Auy vja >o, uu<y\A*vy vfl >o 

for all U 6 15(71), V e D (A*) , and M is a material law in the sense of [3, 6]. More 
specifically we assume that M is of the form 

M (z) = M + zM x + z 2 M {2) (z) 

where is an analytic and bounded L (if, if)-valued function in a ball Be (r, r) for 

some r 6 R>o and Mo is a continuous, selfadjoint and non-negative operator in H. The 
operator Mi 6 L (H, H) is such that 

uM + V\z Mj > c > (2) 
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2 Sturm-Liouville Evolutionary Problems 



for all sufficiently large v e R>o- The operator M (d x ) is then to be understood in the 
sense of the operator-valued function calculus associated with the selfadjoint operator 
3 m (d ) = ±(d -d* ). 

The appropriate setting turns out to be the Sobolev chain 

(H v>k (R, H)) kseZ . 
From [3, 6] we paraphrase the following solution result. 

Theorem 1.3. For J e H vk (R, H) the problem (1) has a unique solution U e H v ^ (R, H). 

Moreover, 

F^(d M (d^+Ay'F 

is a linear mapping in L (H uk (R, H) , H Vjk (R, H)) , k e Z. These mappings are causal in 
the sense that if F e H uk (R, H) vanishes on the time interval ] — go, a], then so does 
(d M (d^) + A)' 1 F, a e R, k e Z. 

Remark 1.4. Note that if U e (R, H) and J e H U)k (R, H) equation (1) actually makes 
sense in H vk -\ (R, H) . Initially the solution theory is for the closure of (doM (d^ 1 ) + A) 
as a closed operator in H vk _\ (R, H) , but it is 

(d M (cX 1 ) +A)U= (d M (do 1 ) +A)U 
with equality holding in H v k _i (R, H) . Indeed, for <\> e H u k (R, H) n D(A*) we have 

(<P\(d M (do 1 ) + A) U)^ = {(dlM (do 1 )' + A*) m\ k _ 1>0 

= <0|M (^) doUy uk _ lj0 + (A*<f>\U\, k -i, 

and we read off that U e D (A) if A is considered in H u>k -i (R, H) (rather than H Vik (R, H)) 
giving 

AU = (d M (d^ 1 ) +A)U- M (do 1 ) d U, 
= (d M (do 1 ) +A)U- d M (do 1 ) U. 

The rigorous argument is somewhat more involved, see [5, 4]. This observation, however, 
motivates dropping the closure bar throughout. 

2 An Application: An Evolutionary Problem Involving a 
Sturm-Liouville Type Operator with an Impedance 
Type Boundary Condition 

We exemplify the outlined theory by a (1 + 1) -dimensional example. 
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Consider 2 

A = 

with an impedance type boundary condition implemented in the domain of A given by 
w J £ H ufi (H, ff (5, /) L 2 (/) © i? (d, I)) | a (d^ 1 ) s - d^v e H ufi (r, i? (3, / 

where H (d, l\ denotes the completion of the space of smooth function with compact 

support in I =] — 1/2, l/2[ with respect to the graph norm of the derivative operator d. 
The space H (d, I) is the domain of the adjoint of 8 also equipped with the corresponding 
graph norm. We focus here on the finite interval case. It should be noted, however, that 
the same reasoning would likewise work for the half infinite interval case, say I = M >0 . 
Indeed this case would be in a sense simpler since only one boundary point would need to 
be considered. We assume 

to be uniformly continuous bounded-analytic-function- valued mapping. As a matter of sim- 
plification we shall assume that {z >— > a (re, z)) xeI is a uniformly bounded family of bounded 
functions, which are analytic in a ball Be (0, 2r) centred at with radius r e R>o. Then 
surely (z >—> a (x, z)) xeI is also a bounded family of analytic functions in Be (r, r) and for 
u > 2r we have a continuous linear and causal mapping: 

a (d^) : H Vt0 (R, L 2 (/)) - (R, L 2 (/)) 

^(t^(i^a(i, ^o -1 ) (t, x))) . 

Assuming further that the distributional derivative a' ■= (x i— > a (rr, ■ )) is such that 

^ a' (re, 2)) xeAJV 

is also a uniformly bounded family of bounded functions, which are analytic in Be (0, 2r), 
N a Lebesgue null set, we get a bounded linear and causal mapping 

a' (d^) : H VtQ (R, L 2 (/)) - (R, L 2 (/)) 

and the product rule 

3 (a (d^ 1 ) s) = a' (do 1 ) s + a (d^ 1 ) ds 
holds for s e D (d). As another simplification we assume that a is real in the sense that 

a (x, z)* = a (re, z*) 



2 This 3 x 3-system has been chosen, rather than an alternative 2 x 2-system formulation, since it shows 
conservativity in a more obvious way, see Footnote 3. 
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2 Sturm-Liouville Evolutionary Problems 



for x e I and z e Bq (r, r). 

Such an operator A combined with a suitable material law yields an evolutionary problem 
of the form 

(d M (do" 1 ) +A)U = J. (3) 
Indeed, we consider a material law operators of the form 3 

//c 
M (do" 1 ) = «! -u^d, 1 

\ /io^o 1 z + V^o 1 + /^o 2 

where e : L 2 (I) — ► L 2 (I), k : L 2 (I) — > L 2 (/), K\ : L 2 (I) — > L 2 (J) are suitable continuous, 
selfadjoint, non-negative mappings and and fj,Q : L 2 (I) — > L 2 (I), \i\ : L 2 (J) L 2 (J), 
77 : L 2 (/) — > L 2 (J) are continuous and linear. We assume of course that the coefficient 
operators are such that (2) is satisfied, i.e. 

ue + £Ke 77 ^ c > 

for some cq g R and all sufficiently large v e !R>o- Such material laws are suggested by 
models of linear acoustics, see e.g. [2], or by the so-called Maxwell-Cattaneo-Vernotte 
law [1, 3] describing heat propagation. In the 1-dimensional case, focused on here, this 
special material law operator can be reduced to a wave or heat equation type partial 
differential operator with a Sturm-Liouville type operator as spatial part. Indeed, assuming 



3 This implies a polynomial type material law operator of the form 
M(d t 



X 1 ) 





-v% 




We see that for kq, ki, e strictly positive continuous linear operators, fi± = and 77 = and A skew- 



selfadjoint, e.g. if a ((9 x ) = 0, we have a conservative system since then M ^ 



and 



Mi 




-/i* 
Mo V 



is skew-selfadjoint making A + Mi and so also yj M 1 (A + Mi) y A/ 1 skew-selfadjoint. Consequently, 



AI Q 1 (A + All ) y M 1 generates a unitary 1-parameter group and "energy" conservation holds in the 
sense that for the solution U of a pure initial value problem we have for t e K>o 



r M Q U (t) 



H 



/ Mo?7(0+) 



H 



or if one prefers to underscore the "energy" metaphor 



E(t) 



<MqU (*) 



H 



M U(0+) =E(0+) 



H 
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additionally that Kq and K\ are strictly positive, two elementary row operations 4 applied 
to 

Ko<5o d 
ki<9 -nl 
d no ed + n + fiido 1 

yield formally 

kq3q d 
Ktdo -/xj 
d„ 1 (edg + r]d + (/i «rVo + A*i) - c^o 1<5 ) 

Clearly, applying the Fourier-Laplace transform to {ed^ + r/do + (/io^rVo + fiA — 8kq d) 
we obtain point-wise, writing VA instead of (im + v), 



(e\ + + (/i «i Vo + A*i) - 8k 1 d^j , 
which for vanishing "damping" rj is indeed a Sturm-Liouville operator 5 : 

rA + q — dpd 



4 



A more common point of view for this operation would be to think of new unknowns being introduced. 

Indeed, if the system unknowns are w then letting 

\v J 

y ■= d^v 

we would get from a line by line inspection of the system 



that 



and so that 



' K d d \ ( s\ / 

Kido -Mo \\w = 



s = - Kq 1 dy 

w = K^ixtv 



ed\y + rjdoy + my + (i a n 1 Vo2/ - dn l dy = f. 



5 If e = and rj has a strictly positive definite symmetric part, i.e. the selfadjoint 9\e rj is strictly positive, 
we arrive at the parabolic type operator 

rjdo + q — dpd 

and writing A instead of (im + v) we get again a Sturm-Liouville type operator 

■q\ + q — dpd, 

where now rj plays the role of r. 
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2 Sturm-Liouville Evolutionary Problems 



with 



Mo^i Vo + Mi 



For our purposes we may allow for general material laws in the problem (3). 

Denoting the inner product and norm of H Ui q (R, L 2 (I) L 2 (I) © L 2 (I)) by ( • | • ) I/ and 
I ' Loo ' respectively, we calculate 




1^,0,0 



v,0,0 



= ^ e ( m o) s|3 (v - d a (d^ 1 ) s )) u0Q + 
+ ^ e <X]-oo,o] ( m o) s|dd a (d^ 1 ) s) uQ o + JKe^al*].,^ (mo) ^ 

= - 9te (^Ix^,^ (m ) (« - 5 a (V) s)X AQ + <He (x^a, (™o) s|dd a (d _1 ) s)^ + 
+ d\t(ds\ X] _ OOt0] {mo)v\ fifl 

= <Re K5x ] _ 00i0] (mo) s|d a (V) s >, 00 + <X]-co,o] ( m o) s|dd a (r^ 1 ) s) u( 
= ^(Xj-oco] K) « ( ■ , +1/2) |5 a (+1/2, c^ 1 ) s ( • , +1/2)\ Q + 
*He<x ] _ 0ifl (m ) s ( ■ , -1/2) |5 a (-1/2, rV 1 ) a ( • , -l/2)\' ■ 



',0,0 



(4) 
(5) 



For this to be non-negative we assume that a is such that 

^Woo.o] (mo) <p\ + c%a (+1/2, <p\ fi > 
for every ip e iJ^i (R). Due to its analyticity at the operators a (x, d^ 1 ) are of the form 

a (d l ) =a + ai d l + a 2 d 2 + d 3 a {3) (d l ) , 
where (Oq 1 ) is bounded. With this we can analyse (5) further. It is 
^ (Xj.co.oi (™o) <p\ + 5 a (+1/2) ip\ fl = 

r° 

= + y> (t)* (d a (+1/2) y>) (i) exp (-2ut) dt 

J— co 



±v | a (+1/2) |</> (t)| 2 exp {-2vt) dt + -a (+1/2) \<p (0)| 2 

-00 ^ 



which is non-negative if we assume 



+a (+l/2) > 0. 



(6) 



12 



Similarly 

^<X]-oo, 0] (™o) <p\±a 1 (+1/2) ip) vfi = 

r° 

= + (£) ai (+1/2) y> (t) exp (-2i/t) dt. 

J— oo 

Assuming that 6 

+va (+1/2) + ai (+1/2) ^ c > 
for z/ g ]R >0 sufficiently large we obtain 

**<*]_«,,«, ( m o) *>| + 5 a (+1/2, d, 1 ) <p\ fi > 
> (±ua (+1/2) + a x (+1/2)) (m ) ^ + 

+ «e ( X] ^ 0] (m ) y>| + ^a® (+1/2, ^) <^ o . 

Due to causality we have 

Ke(x^, 0] (mo) <p\ + <W 2 ) (+1/2, d, 1 ) tp}^ 
= ^^-oco] ( m o) ¥>| + a« (+1/2, rV 1 ) X^, 0] (™o) O ~" V>„ >0 





X] 


-co,0] 


(mo) (/? 


i/,0 


X] 


-oo,0] 


(mo) ^VV 


i/,0 




X] 


-co,0] 


(mo) 


i/,0 


X] 


-oo,0] 


(mo) d 1 x ] 


-oo.o] (™oM,,o 




X] 


-oo,0] 


(mo) (/? 


!/,0 


4 


V 


oo.O] ( m 0) V 


i/,0 



Under this assumption we have 

for all f/ e D (A) if i/ 6 R>o is sufficiently large. Note that by the time-translation 
invariance this is the same as saying 

^{x^ a] {m,)U\AU) ufi ^Q (8) 

for all U 6 -D (A) and all a e R. Letting a — > oo we obtain from this 

JHe^A^^O (9) 



6 If a {d x ) = ao + ai(? 1 it is sufficient to require 

±va Q (+1/2) + ai (+1/2) > 

for all sufficiently large v e K>o- 
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2 Sturm-Liouville Evolutionary Problems 



for all UeD(A). 

We need to find the adjoint of A. It must satisfy 



d" 

| c A* 

dOO 



d' 

5 



in the sense of extensions. We now show that D (A*) is given by 

H (3, 1) © L 2 (I) © H (8, 1)) a ({d i y) s + (do" 1 )* v e H ufi (R, H (d, I 



v 



Indeed, for 



we have 



1 

i\ 3-1 



-a (d l ) d 
Direct computation gives 



o d 

d0 J \-aid, 1 ) d, 



1 

1\ 3-1 

o 



w ]eD(A) 

v 



eH l/fi (R,H(d,I)®H(d,l)). 



d 



1 \ fO d\ ( 

i + 



d x \ ( 8 



o iy \d o 

a (d^ 1 ) 5 


^ -a (d^) \ [08 

o i Uo 



a' (d, 1 ) 




+ 



a'(^)0" 

o or 



Thus, we have 
d 



d~ l 
80 l°° 



+ 



1 a [d^ 


(la{8^) 
\0 

1 a (d, 1 ) 
rV 1 



d 
d 



-a (d^ 1 ) d? 1 



a 1 (ft 1 ) 



d 
d 





o 

0) [v 
1 

a (do 1 j uo 



+ 



dr 1 



+ 



+ 



a' (do 1 ) 
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This implies that 



which is the above characterization of D (A*). Moreover, 

{do){a((d r)(d r)\v 2 )-{ o oju; + 

+ ((V«r n )(S0(;)- 

As a consequence of the similarity between A and ^4* we find by analogous reasoning that 
we have not only (9) but also, indeed more straight-forwardly, 

*fc<y\A*V) vjOJO >o (io) 

for all V 6 D (A*). The calculation is similar to (4) but without the cut-off with ^ (m ). 
Thus we have indeed that 

d M (cX 1 ) + A 

is continuously invertible with causal inverse [p^M (^cT 1 ) + A) 1 : H u j, (R, i7) — ► f/^ (R, i7) 
for every fce Z and any sufficiently large 6 R>o- We summarize our findings in the fol- 
lowing theorem. 
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References 



Theorem 2.1. Under assumptions (2) and (5) we have that for every J e H u k (R ; H) the 
problem (3) has a unique solution U e (R, ff) . The solution operator [d^M (d^ 1 ) + A) 
H Vj k (R, -H") — >• ^,fc (R, -H") continuous and causal for every k e Z and any sufficiently 
large v e M>o- 

Proof. Under the stated constraints the assumptions of Theorem 1.3 are satisfied and the 
result follows. □ 
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